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Abstract 

We study the azimuthal structure of the stellar disks of 18 face-on spiral galaxies, using K'(2.2 / um)- 
band photometry to trace the stellar surface mass density. Assuming the disks are co-planar, we 
characterize their deviation from axisymmetry by the fractional amplitudes, A m (R)/ Ao(R), and 
phases, <p m (R), of the m th azimuthal Fourier components at radii R about the photometric galaxy 
center. We find that most disks exhibit a wealth of non-axisymmetric structures, specifically: (1) 
that about one third of them are substantially lopsided (Ai/Ao 2; 0.20) at 2.5 disk exponential 
scale length, (2) that almost one half of them have strong two-armed spirals with an arm/interarm 
surface-brightness contrasts of order unity, and (3) that typical disks have some intrinsic elliptic- 
ity. We estimate that in the disk plane the characteristic ellipticity of the underlying potential is 
0.045^0 02- However, the spiral pattern couples significantly to the estimate of the intrinsic ellip- 
ticity, and our measurement may represent an upper limit on the "true" potential triaxiality. We 
estimate the radial streaming motions of the disk stars, vr, which are produced by these distortions. 
By averaging over our sample of galaxies and all azimuthal angles, we find (vr) ~ 7 km s _1 due to 
lopsided distortions and (vr) ~ 6 km s _1 due to intrinsic ellipticity. These non-circular motions are 
expected to contribute ~ 0.15mags scatter to measurements of the Tully-Fisher relation. 

Subject Headings: Galaxies: Kinematics and Dynamics, Photometry, Spiral, Structure - Infrared: 
Galaxies 

1 Introduction 

Stellar disks are, by definition, highly flattened. Most galaxy models also assume that their mass 
distribution in the disk plane is axisymmetric. However, deviations from axisymmetry may create 
radial streaming motions and so strongly alter the dynamics and evolution of galaxies. Despite 
their importance, many basic questions regarding the non-axisymmetric nature of disks remain yet 
to be answered. Do disks have a significant intrinsic ellipticity 1 ? Are spiral arms small (<§; 1) or 
large (;> 1) fractional density enhancements? Is the mass distribution symmetric? Recent advances 
in image processing techniques and the development of large-format IR array detectors provide an 
opportunity to re-examine these issues. In this paper, we use these tools to work toward a more 
complete understanding of galaxy disks by measuring the nonaxisymmetric components of the stellar 
mass distribution for a sample of nearby spiral galaxies. 

Much of the previous exploration of the dynamical effects of nonaxisymmetric disks has focused 
on our Galaxy. Various investigators who attribute the observed HI kinematics within the Milky 
Way to a radial motion of the LSR (Blitz and Spergel 1990; Kuijken 1991) conclude that the 
Galaxy's disk cannot be stationary and axisymmetric. Kuijken and Tremaine (1994) arrive at a 
similar conclusion, but without invoking a radial motion for the LSR. Nevertheless, it is difficult to 
determine the precise form of the required perturbation; the resolution of this problem is complicated 
by our position within the very disk whose shape is being analyzed. 

Realizing that it is often easier to solve one's neighbors' problems than one's own, we turn our 
attention to the shapes of stellar disks in a sample of nearby galaxies. Several authors (e.g.Binney 
and de Vaucouleurs 1981; Grosb0l 1985) have used the frequency distribution of isophote shapes in 
order to measure disk ellipticity in external galaxies. These studies rely on the assumption that the 
sample is unbiased with respect to galaxy inclination (an assumption which is probably violated, e.g. 
Huizinga and van Albada 1992), and on the measurement of a single optical wavelength isophote 

1 We use the term "elliptical" throughout, rather than the more precise "triaxial" (a 7^ b and a, b ^> c). Whenever 
we use the term elliptical we mean the intrinsic disk ellipticity (within the disk plane) rather than the observed 
ellipticity arising from projection effects. 
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shape. The distorting effects of dust and star formation on ellipticity measurements are exacerbated 
when only one single isophote is used. 

Our approach has two key advantages over the previous studies. First, we minimize projection 
effects on the apparent disk shapes by selecting objects known a prion to be nearly face-on on the 
basis of their kinematics. Second, we examine the surface brightness distribution of these nearby 
galaxies in the near infrared (K-band, 2.2/rni). It is likely that at this wavelength the surface bright- 
ness distribution provides a direct measure of the nonaxisymmetric features in the mass distribution. 
In the grand-design spiral M51 the K-band light is attenuated by £ 10% even in the most promi- 
nent dust lanes, compared to ~ 50% at I(0.8/um), (Rix and Rieke, 1993, hereafter RR93). RR93 
also found that despite vigorous star formation in the arms, the luminosity of young supergiants 
comprises a small fraction of the total K-band light in M51. Because longer wavelength images 
start being contaminated by emission from hot dust, K-band imaging of face-on galaxies is the most 
straightforward way of assessing the stellar surface mass density of galaxy disks. 

We concentrate on three specific distortions: elliptical disks, lopsided disks, and two-armed 
spirals. The ellipticity of disks may bear the imprint of details of the galaxy formation process. 
Simulations suggest that the galaxy formation process is likely to lead to significantly triaxial halos 
(e.g. Warren et al. 1992), which may harbor elliptical disks. At radii where the halo dominates 
the potential, the shape of the disk is straightforwardly related to the shape of the halo potential. 
Furthermore, the shape of the disk affects the Tully-Fisher (Tully and Fisher 1977) relationship 
between a galaxy's luminosity and HI linewidth (hereafter referred to as the TF relationship). This 
relationship is currently a key tool in the study of the local Universe and its small "intrinsic" 
scatter attests to the homogeneity of the galaxy formation process. Franx and de Zeeuw (1992) 
demonstrated that even a relatively small amount of disk ellipticity (e £ 0.1) introduces significant 
scatter into the TF relation and argued in turn that the observed TF scatter provides an upper 
limit on the ellipticity of the HI disks. The second type of distortion, lopsidedness, may be a clue 
to recent infall activity or interactions and may affect disk evolution by generating large radial 
motions. Finally, measuring the mass amplitude of the spiral arms provides a direct indication of 
the driving force affecting various processes associated with spiral arms, such as gas compression 
and star formation. The fractional amplitude of the driving mass perturbation is generally assumed 
to be less than unity (see e.g. Athanassoula 1984, and references therein), even though there is 
little support for this assumption from the observations (e.g. Schweizer, 1976; Jensen, Talbot and 
Dufour, 1981; Elmegreen and Elmegreen, 1984, 1985; Elmegreen et a/.1992, Elmegreen et a/.1993; 
Grosb0l, 1993; RR93). However, as of now there is no definitive measurement of the dynamical 
strength of spiral arms because the presence of dust and star formation in the spiral arms complicate 
this measurement. A detailed IR surface photometry study of M51 (RR93, see also Schweizer 1976, 
Elmegreen and Elmegreen, 1984) confirmed that the spiral arm mass perturbation is large (fractional 
density enhancements of a factor of 2 to 3) in that galaxy. While that result cannot be generalized 
because of the strong interaction between M51 and NGC 5195, it may indicate that spiral arm mass 
amplitudes have been significantly underestimated. 

This paper is organized as follows. In Section II we discuss the galaxy sample, the observations, 
and the data reduction. In Section III we describe the azimuthal Fourier decomposition of the galaxy 
images, present an overview of the empirical results, and describe the model fitting. In Section IV we 
also discuss the connection between surface brightness and surface mass, the connection between the 
isophote shapes and the potential distortion, ellipticity of spiral galaxy disks, lopsidedness in galaxies, 
and the strength of spiral arms. Throughout Section IV we discuss briefly the implications of the 
observed asymmetries on disk dynamics. Section V presents the conclusions. Appendices provide 
details on the model fits, the derivation of our estimated distribution of intrinsic disk ellipticities, 
and a discussion of the relative contribution of young stars to the total light in the near-IR. 
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2 The Data 



2.1 Sample Selection 

The measurement of disk shapes is simplified by obtaining a sample of face-on (zero inclination or 
i = 0°) galaxies. Inclinations are generally determined from a galaxy's apparent axis ratio under the 
assumption that disks are intrinsically round. To measure intrinsic disk shapes we obviously must 
use a shape-independent means of estimating the disk orientation. We attempt a purely kinematic 
determination of the inclination, i, by "inverting" the TF relation. In its standard form, the TF 
relation is given by 



with = -19.55 mag, a = -7.48, and v* = 158 km s" 1 (Pierce and Tully 1992), where Wl, 
represents the width of the HI emission profile at the 20% of the maximum level, corrected for 
the internal velocity dispersion of the gas. With this relation we assign to each galaxy of absolute 
magnitude Mb a " characteristic velocity", v(Mb) (which approximates the disk circular velocity). 
Inverting Eq. (1) to solve for i results in 



Here we determine a galaxy's inclination by comparing its observed HI line width with its char- 
acteristic velocity, derived from its absolute magnitude (This magnitude, in turn, is calculated by 
assuming pure Hubble flow and Ho = 75 km s _1 Mpc -1 .) Note that for nearly face-on inclinations 
the observed linewidth is very sensitive to i, while the projected shape is not. For example, if a 
face-on galaxy with i> c ; rc = 200 km s _1 is tilted by 15°, its observable linewidth has already increased 
to 104 km s _1 , about a quarter of the true value, while the axis ratio of the projected circular disk 
has changed only to 0.97. However, even galaxies with perfectly face-on stellar disks may have a 
linewidth larger than the vertical gas velocity dispersion (Lewis 1984), for example because many 
gas disks are warped at large radii. As a consequence luminous galaxies with narrow observed HI 
linewidths are much rarer than expected from a uniform distribution in cos i and an equation anal- 
ogous to Eq. 2 (Lewis 1984). We did not restrict ourselves to well known examples or samples of 
face-on galaxies (e.g. Lewis 1984) because many of these studies are based on a stringent axis-ratio 
pre-selection of the objects. Instead, we searched the largest databases available, the Huchtmeier- 
Richter catalog (Huchtmeier and Richter 1989) for HI linewidths and the CfA Redshift database 
(Huchra 1993) for apparent galaxy magnitudes and redshifts, in which most entries are unrelated to 
the question at hand (therefore avoiding any morphological bias). For all candidate galaxies with 
published HI profiles, we examined the original data to assure that they were of good quality and 
that there was no evidence for a gas-rich companion in the beam. We examined images on POSS or 
ESO sky survey images to restrict our sample to isolated field spirals with i £ 25° . 

We list in Table 1 our final sample of 18 galaxies. We present the galaxy's name in column (1), 
its type in column (2) from the Revised-Shapley Ames catalog (Sandage and Tammann 1981), its 
apparent blue magnitude and diameter from RC3 (de Vaucouleurs et a/.1991) in columns (3) and 
(4) respectively, its absolute magnitude corresponding to Ho = 75 km s _1 Mpc -1 in column (5), its 
recessional velocity and linewidth from the Huchtmeier-Richter HI catalog corrected for an internal 
velocity dispersion of 10 km s _1 (Tully 1988) in columns (6) and (7) respectively, its inclination 
from Eq. 2 in column (8), and the total K' exposure time in column (9). The median luminosity 
of galaxies in this sample is 1 . li* and their median distance is 31 Mpc. Because the catalogs used 




(1) 




(2) 
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to identify the sample have heterogeneous selection criteria, our resulting sample is not statistically 
complete. Nevertheless, there is no obvious correlation between a galaxy's inclusion in the catalogs 
and the details of its stellar disk shape. In particular, all galaxies with (apparent) axis ratios between 
0.85 and 1 should stand a comparable chance of entering the catalogs. Therefore, we do not believe 
that a bias toward either circular or slightly elliptical disks has been introduced into the sample. 

2.2 Observations and Data Reduction 

We present K'(2.2 / um)-band images for all sample galaxies in Figure 1. These are complemented 
by additional I(0.8//m) images for six sample galaxies (see Table 1), which were used to construct 
external error checks. The K'-band data were obtained during 1992 Nov 11 through 14 at the Las 
Campanas 2.5m telescope with a NICMOS3 256x256 array. The K'(2.2 / um) filter has an effective 
wavelength of 2.16/um and a bandwidth of 0.33/um and is designed to reduce the background sky 
flux compared to the standard K filter by lowering the long wavelength cutoff of the bandpass. 
The image scale is 0.405" pixel -1 with a field-of-view of 1.73' x 1.73'. The galaxies have angular 
sizes comparable to the field-of-view and the effect of this on the sky determination is discussed 
in §3.4.1. The atmospheric transparency was good but not photometric and no flux calibration 
was attempted. During the four nights, the seeing at K' was £ 0.8". After multiple images were 
combined, the effective seeing was typically between 1" and 1.2 " . The total observing time for 
any given galaxy was split into about 10 "on" and 10 "off" exposures, alternating every few minutes 
between the object and a blank patch of sky ^ 5' away. Small random components were added 
to these offsets to improve flat fielding and to avoid the accumulation of signal in images of faint 
objects in the "off" images. The total exposure time on each object is listed in Table 1. Additional 
5 minute I-band exposures were obtained for six objects using the Las Campanas lm telescope and 
a Tektronix 2048 x 2048 CCD. 

The reduction of IR array data differs from the reduction of CCD data principally because the 
background signal at K' is much larger than at optical wavelengths. The sky flux is so overwhelming 
that the detector must be read out at least every 30 seconds even in the K' band. At a radius 
of three disk scale lengths a galaxy disk is already 8 magnitudes below the sky level! Therefore, 
small errors in the sky subtraction can introduce systematic errors and increase the photon noise. 
To achieve the desired level of precision, images of "blank" sky are used for both flat-fielding and 
background subtraction. If, as we did during our observations, equal amounts of observing time 
are spent on object and on sky, then straightforward sky-subtraction will increase the noise by 
\/2. This degradation can be mitigated by using a sky image that is the combination of the "off" 
frames of several objects. Unfortunately, the sky-subtraction noise cannot be reduced arbitrarily 
in this way, because the sky flux at K' is due both to OH emission in the upper atmosphere and 
thermal emission in and around the telescope. These two sources illuminate the detector somewhat 
differently and are time variable. Consequently the total illumination pattern changes shape and 
scale as the temperature and the atmospheric conditions change. This variability limits the total 
time over which a sky signal can be combined usefully to less than two hours. We typically achieve 
a sky flatness good to 0.03% across the frames. Standard IRAF 2 tasks were used to flatten, align, 
combine, and sky subtract the images. 

2 IRAF is distributed by the National Optical Astronomical Observatories, which are operated by AURA, Inc., 
under contract to the NSF. 
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3 Data Analysis 



3.1 Fourier Decompositions 

Fourier analysis is a natural way to analyze the non-axisymmetric component of the light in a face-on 
galaxy. It has the advantage over model fitting (see §3.2) that no prior assumption is made about 
the intrinsic light distribution. Fourier analysis techniques have been used frequently in the past 
(cf. Grosb0l 1987, Elmegreen et a/.1989, Elmegreen et a/.1992, RR93, Elmegreen et a/.1993, and 
references therein) and we summarize only briefly the particulars of our implementation. We fix the 
center of a polar coordinate system at the brightest point of the galaxy in K' , presuming that it 
accurately represents the bottom of the overall potential well. This choice is justified a posteriori, 
because in their inner portions the sample galaxies only show small asymmetries about the chosen 
center. Once the center is set, the surface brightness distribution distribution iJ,(R,ip) is expressed 
as a Fourier series: 

oo 

»(R, <p)/(ii(R)) = E A m (R)e^-^\ (3) 

m = l 

where <p denotes the azimuthal angle, m the azimuthal wavenumber, and A m and <p m are the 
associated Fourier amplitude and phase, respectively. The average surface brightness at radius R is 
given by {n{R)). 

The analysis sequence consists of several steps. First, we mask all detectable point sources. 
Second, we re-bin the image onto a (R, <p) grid, using 30 bins in radius and 24 bins in azimuth. 
The rebinning leads to a large speed-up in all subsequent calculations. The minimum and max- 
imum radii for the analysis were set to ~ 6 pixels and 170 pixels, respectively. At small radii 

lOpix ) some interpolation is unavoidable and the data on the polar grid are not fully inde- 
pendent. At larger radii where the number of original pixels per new bin is large, we use the 
median of the pixels to obtain //(i?, <p). The average and median produce indistinguishable results. 
T he error at each grid point, Ajj,(R,ip), was taken to be the quadrature sum of the Poisson noise, 
^/((source) + (sky))/(number of pixels), and the "flat-fielding error", £fi a t ■ (sky). The flat-fielding 
error also incorporates sky subtraction errors, which chiefly arise from imperfect flat-fielding. Even 
though the NICMOS data can be flatfielded extremely well (f/; a t ~ 3 x 10 -4 ), flat-fielding errors 
become dominant at large distances from the galaxy center. In this regime, the errors in Figures 2 
and 4 become spatially correlated because the flat-fielding error is assumed to be the same across the 
whole frame. Lastly, after a (R,<p) grid has been established for \i and A/j,, the Fourier amplitudes 
and phases are determined by a least squares fit to the one-dimensional array of data at each R. 
Solving the least squares problem by means of Normal Equations (cf. Press et al. 1986) also yields 
the variances for each parameter. Systematic uncertainties in this process are described below where 
we compare the results for an identical analysis of the independent I and K' images (§3.1.4). 

In Figure 2 we plot selected Fourier components for the 18 galaxies. These components were 
chosen (from the complete Fourier expansion) because they are the only ones that differ significantly 
from zero in most sample galaxies. The 0th order Fourier amplitude, Ao, has no phase associated 
with it and merely reflects the mean flux in the annulus, (^(R)}. It is plotted as a function of radius 
in the top panel. For most galaxies the surface brightness profile of the disks is well approximated 
by an exponential (a straight line in these plots) and we have chosen the exponential scale length as 
the radial unit of the plots. Our surface photometry typically extends to about three scale-lengths 
(as determined in Section 3.2 and listed in Table 2). We present an overview of the results for the 
various m > 1 Fourier components in the subsequent sections. 
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3.1.1 m=l 



Panels two and three in each column of Figure 2 show the m = 1 amplitude and phase. The 
amplitude Ai(R), characterizing the fractional variations in the surface brightness at radius R, is 
related to the fractional offset, AR/R, of the isophotes from the center by 



For an exponential disk with scale length R exp this implies A^f/j, = R/R exp ■ AR/R. 

In the inner disk (1 to 2 scale lengths) Ai distortions are generally weak. However, beyond these 
radii a sizable fraction of galaxies exhibit significant lopsidedness, Ai ^ 0.2. The most prominent 
examples are IC 2627, NGC 1309, NGC 1325, NGC 1642, NGC 2485, and NGC 6814. An inspection 
of sky survey plates did not reveal any obvious environmental differences between these lopsided 
galaxies and the others. One notable exception is NGC 1309 which has an apparent companion 
visible in our I-band CCD image. A histogram of the lopsidedness in the sample, measured at 
2.5 disk scale lengths, is shown in Figure 3. Because the amplitude of the "true" lopsidedness, 
Ai(2.5R exp ), is a positive definite quantity, the expectation value of its measurement, A\(2.'oR exp ), 



in the presence of an error AAi is y Af + AAf (see Section 4.2.2 for a more extensive discussion). 

We have accounted for this when plotting the histogram. The sample mean of A\(2.'oR exp ) is 0.14; 
this value drops to 0.11 if the largest value is excluded. We chose R = 2.bR exp as a radius of 
comparison because this is the largest radius for which well determined estimates of Ai exist for 
most the galaxies in the sample. 

It is important to note that these results are not affected by an improper choice for the center of 
our polar coordinate system. If the galaxy isophotes were concentric, but we had chosen the wrong 
center, then Ai(R) would diverge as l/R for R — > 0, as long as there is a finite intensity gradient. 
Furthermore, the asymmetries cannot be a consequence of projection, because that operation is 
bi-symmetric. Finally, for nearly face-on galaxies the dust extinction on the "near" and "far side" 
is indistinguishable and will not induce an asymmetry. 

3.1.2 m = 2 

The fourth and fifth panels from the top in Figure 2 show the m = 2 amplitudes and phases for all 
sample members. The inclination selection of our sample (i ^ 20°) limits the apparent ellipticity due 
to projection effects for circular disks to A'^iR ~ 2R exp ) ^ 0.06. However, the majority of objects 
exhibit y^'s much larger larger than 0.06 at some radii. Such large amplitudes can arise from either 
disk intrinsic ellipticity or two-armed spirals. Elliptical disks (or bars) and spirals arms differ in 
their pitch angle, 6 = atan(91ogi2/9y), which is expected to be ~ 90° for bar-like distortions, and 
0° < 6 < 30° for spiral arms. Therefore, the relationship between <p and R provides information on 
the nature of the structure. Not surprisingly, the strongest m = 2 elliptical distortions are found in 
the central bars of the galaxies. The most notable examples in the sample are NGC 600, NGC 1015 
and NGC 6814, which have large A'^s for R < 50" at constant <p (or dlogR/dcp = oo). At radii 
beyond this, the strongest features are spiral arms, which have amplitudes of 0.15 <> A^ 0.6. In 
many of the galaxies in the sample, the arms wrap continuously through at least 180°; in IC2617, 
NGC 1642, NGC 1703, and NGC 6814 they wrap at least 360°. 

3.1.3 m = 3,4, and 6 

The bottom three panels in each column of Figure 1 contain the amplitudes of the third, fourth, and 
sixth Fourier components. Note that the vertical scales in these panels have been stretched by 5/3 



Afi 



Slog// Ai? 
dlogi? ' ~R' 



(4) 
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relative to the m = 1 and m = 2 panels because the higher frequency amplitudes are significantly 
smaller than Ai and A2. Notable exceptions to this are found if strong and thin bars are present 
(NGC 600 and NGC 1015), then A4 and Aq are comparable to A2. If a bar were infinitely narrow, 
such that I(<p) = + Hv + 180°), then all even numbered coefficients, A2 n (n = 1, 2, ...), would 
be equal. In practice, bars have finite thickness and the higher order Fourier amplitudes fall off. 
In some two-armed spirals the non-axisymmetric power is concentrated predominantly in m = 2, 
(^2} 2i 2(^4), which implies that the spiral arms are quite broad (e.g. ESO-436, NGC 1309, 
NGC 1703, NGC 2466, NGC 2485). However, the two galaxies with the most prominent spiral 
patterns in the sample, IC 2627 and NGC 6814, also show large m = 4 amplitudes, because the 
strong two arm spirals are narrower than a simple cos 2<p variation. In most cases, terms of order 
m > 4 are small, e.g. Aq <C 0.1, for nearly all spirals. Three galaxies in the sample, ESO-436, 
NGC 1376 and NGC 7309 have three-arm spiral patterns that can be easily discerned in the images. 
These patterns are clearly reflected in the significant m = 3 Fourier amplitudes. 

3.1.4 Comparison of I and K' Band 

Because other K' data of this type are unavailable, the most direct external error check we can 
perform is a comparison to our I-band CCD images. However, images in the two colors are not 
expected to be identical. A multi-color study of M 51 (RR93) demonstrated that the K' and I band 
Fourier decompositions can differ, mainly due to obscuration by dust. Nonetheless, many properties 
such as the overall variations of the m = 2 amplitudes and their phases provide a benchmark on the 
reliability of the K'(2.2 / um) data. 

We obtained large field, 20' x20', CCD images in I for six sample members: NGC 1302, NGC 1309, 
NGC 1376, NGC 1642, NGC 1703 and NGC 2466. We transformed these CCD images to the format 
of the K' images and estimated their sky level with the same procedure used for the NICMOS 
data. We then performed a Fourier decomposition using the same grid as that used for K' data. 
The resulting comparison for the six objects is shown in Figure 4. The top panel contains the I-K' 
color, subject to an arbitrary offset. All galaxies become slightly bluer with increasing radius for 
radii <, 40". The interpretation of this color gradient in terms of population and dust extinction 
gradients is beyond the scope of the present paper. At radii ^ 40", the uncertainties in the sky level 
increase beyond ~ 0.1 — 0.2 mag. However, for all galaxies the colors at large radii are consistent 
with an extrapolation of the color gradient at smaller radii. We take this as an indication that we 
have conservatively estimated our sky errors. For most of the galaxies, the m = 1 and 2 amplitudes 
and phases are consistent over most of the radial range. Wherever the Fourier phases differ by large 
amounts, differ either by 360°/m or in regions where the associated Fourier amplitude is very small. 
Note that discrepant amplitudes do not primarily occur in the outer low signal-to-noise portion 
of the profile, but rather in the inner, well measured portion of the profile. These differences are 
presumably caused by dust. Note in particular that the agreement of the m = 1 distortion between 
I and K' shows that these asymmetries are neither due to gradients in the sky background nor due 
to asymmetric dust extinction. For 20" < R < 60", m = 2 is almost always larger in the K' data 
than in the I data. This is consistent with the results for M51 (RR93), where a significant fraction 
of the I-band m = 2 amplitude is suppressed by dust along the ridges of the arms. 

The principal function of the I-band images is to test for distortions in the K' images caused by a 
poor sky level determination whenever a galaxy fills the field-of-view. For this test we compared the 
luminosity profiles derived from the trimmed I images with the results from the untrimmed images, 
for which the sky level was determined from the much larger field-of-view. In all cases the two sets 
of profiles agree within the error bars. Note that this test does not depend on the absolute sky level 
(which is much higher in K' ), but only on the degree to which the galaxies fill the field-of-view. 
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The general agreement between the K' and I results, and between the results from large and 
small field-of-view I images, substantiates the assigned K' error bars. 

3.2 Model Fitting 

Because the m = 2 amplitude from the Fourier analysis reflects a mixture of inclination effects, 
bulge and disk ellipticity, and spiral arms, it cannot be used to place unambiguous constraints 
on the shapes of disks. Without further information these various effects cannot be untangled 
rigorously. We make progress by concentrating on the ellipticity of the disk at a radius-independent, 
but arbitrary, position angle. By fitting a specific model to the images averages radially over the 
distortions created by the spiral arms. As we derive in Section 4.2.1, for radii greater than ~ R eX p the 
isophotes of a disk residing in a non-rotating logarithmic potential of ellipticity e$ (in the disk plane) 
have an ellipticity, e 8S0 (i?) = e$(l + R exp /R) (Franx and de Zeeuw, 1992). The first step in this 
analysis is to derive the radial range over which the disk dominates the total light (/J, disk > ^/J-buige) 
by fitting an axisymmetric bulge-disk model to the whole image. Then we estimate e$ by fitting 
a disk model with an ellipticity profile given by e 8S0 (i?) over a restricted radial range of the image 
(see Appendix A for details). We adopted this procedure, rather that fitting simultaneously non- 
axisymmetric bulges and disks of independent ellipticities and position angles. This simultaneous 
fit often yielded strongly distorted bulges and disks oriented at right angles. The superposition 
of a perpendicularly oriented bulge and disk produces a much less flattened overall image than 
each component separately: a presumably unphysical result. For each galaxy our fit provides an 
estimate of the scale length, R exp , the ellipticity of the potential, e$, and position angle, PA, and 
their associated uncertainties. These results are presented in Table 2. With a few exceptions (e.g. 
NGC 1302 and NGC 2718), the values for are found to be £ 0.1. 

Because a perfectly open (a pitch angle of 90°) spiral arm is indistinguishable from disk ellipticity, 
we use simulated data, constructed from a perfectly axisymmetric and exponential disk with a 
superimposed spiral, to study the coupling between the arms and the disk ellipticity. The simulated 
spiral winds logarithmically with a pitch angle of 20°, and extends from one to three disk scale 
lengths. These parameters are characteristic of the spiral pattern found in our sample. We set 
for the simulated galaxy to match the one observed in NGC 1703 and assumed A^{R) = A2(R)/2. 
The simulated galaxy image is shown in Figure 5, and seems to mimic a typical grand design spiral. 
We then add noise and sky errors to this image that are matched to the noise properties of our 
data. Finally, we determine the disk's ellipticity using the same fitting procedure as before and 
find e = 0.03 for this simulated spiral. This result suggests that although spiral arms do contribute 
significantly to the measured disk ellipticity, they do not dominate the signal in the e$ estimates 
for most sample members. In our final analysis we treat this interplay between disk ellipticity and 
spiral arms by quadratically adding an extra error term of Aesp = 0.03. 

4 Discussion 

4.1 Surface Brightness vs. Surface Mass Density 

To assess the implications of the K'(2.2//m) Fourier decomposition and of the model fitting, it is 
important to understand the degree to which surface brightness traces the surface mass density. 
This mapping is complicated by dust extinction and by spatial variations in the stellar mass to light 
ratio. The first effect is largely eliminated by observing in the K' band, where the dust optical depth 
is reduced by a factor 6 compared to 7(0.8//) (e.g. RR93). But because a population of young stars 
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is brighter at all wavelengths than corresponding old stars, the common lore that "red light arises 
from old stars" may be dangerously misleading. 

There are two ways of estimating the contribution of young stars to the surface brightness en- 
hancement in the spiral arms. First, because young stellar populations are both brighter and bluer 
than older ones, one can compare the arm/inter-arm color change to the arm/inter-arm brightness 
contrast (Schweizer, 1976; Jensen, Talbot and Dufour, 1981; Elmegreen and Elmegreen, 1984; Ken- 
nicutt and Edgar, 1986). These authors demonstrate that the color change across the spiral arm is 
too small to explain the brightness contrast merely through azimuthal mass-to-light ratio variations; 
consequently they argue for a substantial enhancement of the stellar surface mass density in the 
spiral arms of grand design spirals. Alternatively, one can estimate the contribution of luminous 
young stars in the IR directly. RR93 attempted to do so by measuring the CO band-head strength 
in the disk of M51 and concluded that although young stars may contribute significantly (~ 20%) 
to the K-band light in the arm crests, they cannot be responsible for the bulk of the light from the 
spirals arms. The situation for M 33 appears similar. Regan and Wilson (1993) obtained B and V 
band photometry for part of the southern spiral arm of M33 and detected red supergiants in some of 
the stellar associations. However, even in those regions, the RSGs comprise only 10% of the K band 
flux (assuming that they are 2 magnitudes redder in V-K than the underlying disk). Therefore, at 
least where one can study the spiral arm population in detail, the light from young stars contributes 
only between 10 and 20% of the total K-band light. 

Because the galaxies in our sample are too distant for a detailed study of their stellar content, 
we supplement the previous discussion by a third, independent argument: a brief discussion of the 
theoretically expected "age distribution" of the K-band light. As detailed in Appendix C, we use 
the stellar population models of Chariot and Bruzual (1991) and simple assumptions about the star 
formation histories of galaxies (e.g. Kennicutt 1983; Kennicutt, Tamblyn, and Congdon 1994) to 
estimate the fraction of the total K light contributed by stars younger than age t y , Lk{< ty) / 1 Lk{< 
to), where to is the current age of the Universe. This age distribution is shown in Figure 6 and 
it implies that 8% of the IR light arises from stars with t ^ t 0T ut ~ 10 8 yrs and 18% from stars 
with t I0t or i,n. Because stars phase mix or drift out of the spiral pattern in ^ I0t or i,a, these 
models imply that the upper limit on the total nonaxisymmetric power (all Fourier components 
except Ao summed) that could arise from azimuthal population variations is 20%. However, this 
line of reasoning also suggests that about half of the light from stars with t < I0t or i,a should come 
from stars with t < t or i,a. Because stars with ages substantially less than t or i,a have not drifted 
far from their birth place, the azimuthal population variations arising from these stars should be 
sharply defined. When the distribution is sharply defined in azimuth (as in the bar of NGC 1015) 
then A2 = A4 = Ag. If the spread in azimuth is less than 27r/m n , we expect this equality to hold 
for all even A m , m <C m n . However, A4 is observed to be generally considerably smaller than 
A2 (see §3.1.3), in conflict with these model predictions. The discrepancy may be removed if the 
models are overestimating the K' luminosity of very young stars. Nevertheless, the results from these 
calculations suggest an upper limit to the contamination of the total IR luminosity from young stars 
of 20%. 

The empirical evidence and the model predictions suggest that the assumption of mass tracing 
light at K' is valid within ~ 20%. The regime in which it is most likely to fail is at large azimuthal 
wave numbers (m > 2) and small amplitudes (A m < 0.1). 

4.2 How Elliptical Are Galaxy Disks? 

Are galaxy disks triaxial? Or equivalently, are they elliptical when viewed perpendicular to the 
disk plane? We cannot answer these questions as straightforwardly as we can measure lopsidedness, 
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because two-armed spirals and projection effects affect measurements of triaxiality. Therefore, we 
focus on a specific form of disk ellipticity: a disk residing in a logarithmic (halo) potential of constant 
fractional distortion and without figure rotation. The choice of a nonaxisymmetric and logarithmic 
potential is motivated by the approximately flat observed rotation curves of spiral galaxies out to 
several optical radii (e.g. Begeman 1987 and references therein). Setting the pattern speed to 
zero is motivated by Occam's razor and by the fact that it avoids dealing with resonances. These 
assumptions do not necessarily conflict with the existence of wound spirals arms of finite pattern 
speed, which may occur in a triaxial potential. 



4.2.1 Isophote Shapes and the Gravitational Potential 

With the assumption that most stars at R ^ R eX p are on loop orbits one can derive a unique relation 
between the shape of the potential and the shape of face-on isophotes (see Franx and de Zeeuw 
1992). Starting with a flattened potential with (small) ellipticity e$ in the equatorial plane and a 
characteristic circular velocity v c , 

2 2 

<I>(R,f) = Y lnR2 ~ yy cos + const., (5) 

and applying first order perturbation theory (as outlined in Binney and Tremaine 1987, pg. 146) 
one can obtain the set of closed loop orbits at the mean radius i?o 

i?(^) = i?o(l-ycos2^), (6) 

v v (f) = v c (l + £$ COS lip), (7) 

and 

vr{<p) = -« c e$ sm2f. (8) 
These orbits have the same ellipticity as the potential, but are elongated (f^tj = 90°) perpendic- 
ularly to the major axis of the potential (fmaj = 0°). The velocity along the orbit varies by twice 
e$, with Vp^max) occurring at R(min), while the extrema in vr are off-set by 45°. An ensemble of 
stars on such orbits must satisfy the continuity equation 

^[R-»(R,<p)-v R (<p)]+-^[rtR,<p)-v ip (<p)\ =0, (9) 

where \i denotes the surface mass density. If the radial surface brightness profile is an exponential, 
we can write \i(R, f) as 



H(R, if) = Ho exp 



R £iso(R) „ \ 
(1 cos 2f) 



Rexp 



(10) 



where /io is the central surface brightness and e 8S0 is the ellipticity of the isophote. We obtain from 
Eqs. (6) - (10) 

ei.o(R) = e*(l + ^f^. (11) 

Eq. (11) demonstrates that the disk isophotes are always more distorted than the potential, e.g. by 
a factor of 1.5 at two scale lengths. The divergence of e 8S0 (i?) as R approaches zero is of no practical 
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consequence because the assumption that stars are on loop orbits fails for R <C Rexp (e.g. Kuijken 
1993); we only consider the regime for which R ^ R eX p- 

So far we have only considered streamlines of closed orbits, such as those of a kinematic tracer 
with negligible radial velocity dispersion (e.g. H I). Stars, however, have a non-negligible dispersion, 
and move on orbits that are "sired" by the closed parent orbit. The general form of a non-closed 
orbit in the above potential is 

R(t) = R (l - y cos 2SM + C* cos V2SM) (12) 

(Binney and Tremaine 1987), where £1 is the angular frequency of the unperturbed orbit. As evident 
from Equation 12, the stars oscillate with an irrational frequency about the location of the parent 
orbit. To first order, the amplitude C of this oscillation, which can be derived by differentiating Eq. 
12, is given by C 2 = jv 2 c . In the solar neighborhood C ~ 0.2 and hence is much larger than the 
expected contribution from e$/2. An ensemble of such orbits forms an outer envelope of ellipticity 
tauter = 1 - (1 - ^ + C*)/(l + ^ + C*) and an inner envelope of e inner = 1 - (1 - ^ - C)/(l + $f - C). 
Accounting for the greater weight of the outer parts to the observed isophote (due to the slower 
azimuthal velocities at apocenter) and averaging the inner and outer edge of such an annulus, one 
finds that the ellipticity of the isophotes only changes by — e$C 2 compared to the closed orbits, 
with all lower order terms vanishing. Therefore, the isophote shapes are independent of the stellar 
velocity dispersion in the limit uj; < u c . We will use this approximation throughout the paper and 
it enables us to constrain the potential shape simply from the isophote shapes. 

4.2.2 Estimating the Distribution of 

In Section 3.2 we described how we estimate the observed disk ellipticity e$ and its uncertainties 
for each of the sample members. If our galaxies were perfectly face-on and the ellipticity error 
were negligible, the results from Section 3.2 would give an immediate observational estimate of the 
intrinsic ellipticity distribution, P(e^), for the sample. 

In practice our estimate of must account for two complicating effects: ambiguities in the 

deprojection and errors in ellipticity measurement. When the disk has a small, but finite, inclination, 
there is a range of intrinsic ellipticities that can result in the observed ellipticity. If we denote 
the projected position angle of the intrinsic disk major axis as PA ma j and the projected position 
angle of the disk normal direction as PA n , then this projected ellipticity depends on PA ma j- n = 
PA ma j — PA^. Because our inclination estimates come from single beam H I measurements, we 
have no prior knowledge of PAf, we therefore must assume that PA ma j- n is uniformly distributed. 
The ellipticity errors, Ae enter the estimate of -P(e$) in two ways: (1) can obviously not be 
measured perfectly, and (2) measurement errors lead to a biased estimate of because ellipticities 
are positive definite quantities. Fortunately, analogous probability distributions for positive definite 
quantities have been derived in other astrophysical contexts (e.g. Wardle and Kronenberg 1977). 
We describe the application of those distributions in detail in Appendix B. 

In the subsequent sections we describe how to estimate from our sample in both a para- 

metric and a non-parametric fashion. In the first approach we do not need to adopt a functional 
form for P(e^). This is an important asset because we have no physical arguments for any partic- 
ular form. However, such a technique will not recover the "true" distribution even for an infinite 
number of measurements. In contrast, we can obtain an unbiased estimate if we adopt a parametric 
functional form for the P(e^). Our philosophy here is to first estimate non-parametrically 
and then to use this result to choose a functional form for the parametric estimate. 

a) Non-Parametric Estimates 
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As shown in Appendix B, one can determine from the observed parameters for each galaxy j, 
the probability distribution of intrinsic ellipticities, Pj(e|,), cost, e$, Ae^. This estimate is based 
on the assumption that the angle specifying the disk orientation at given inclination is distributed 
uniformly. The probability distribution for the entire sample is then 

p(4) = j^ E (13) 

j = l,N 

This distribution is shown as the thick solid line in Figure 7. We estimate the 95% confidence region 
for this distribution due to our finite sample size by means of "bootstrapping" (e.g. Press et at, see 
Appendix B). These estimates are shown as the two dashed lines enveloping our estimate of P(e^). 

This estimate will never converge to the true answer. Consider for example an ensemble of 
axisymmetric disks, Ptme(e$) = ^( e $); seen at small but finite inclinations. Then all disks will 
project into slightly elliptical shapes. If we have no prior knowledge about Pt™e(f$), we can never 
exclude the possibility that some disks are truly elliptical. Finite measurement errors lead to a similar 
effect. Therefore, our derived distribution P(e$) should be compared to the distribution, Po(e^), 
that would result if we observed a set of perfectly round disks with inclinations and uncertainties 
corresponding to the disks in the sample. The distribution Po and its 95% confidence regions are 
shown as the dotted lines in Figure 7. Po can be thought of as a "point-spread function" because it 
reflects the estimate of P(e$) for an input (^-function. 

The observed P(e$) appears broader than Po(e^), even when considering the range of the con- 
fidence limits. The median for axisymmetric disks is 0.025 but is 0.07 for our sample galaxies, 
showing that the underlying stellar disks in these galaxies are slightly elliptical. However, these 
non-parametric estimates do not provide a good measure for a "characteristic" ellipticity, to, and 
complicate the presentation of confidence limits. 

b) Parametric Likelihood Estimates 

The observed P(e^), displayed in Figure 7, prompted us to chose an exponential as the parametric 
trial function for P(e^): 



P(e$, eo) = — exp 

to 



(14) 



Our task is now to determine the most likely value of to • Note that the mean and the variance of this 
distribution are both equal to cq. For each assumed e^, inclination, and PA ma j- n , we can determine 
the expected ellipticity for galaxy j in the noiseless C&S6, ^no — noise- 

In the presence of measurement 
error, Ae, the distribution of measured e$ for galaxy j is given by (Wardle and Kronenberg 1977) 

n / \l e $ t I t^^no — noise \ / ^ <3> ' ^no — noise \ / -i r \ 

Pj(^)\ei,cos 3 (i),PA m ^ A = ^ h I ^ j exp I — I , (15) 

where 7o is the modified Bessel function. Integrating Pj over all (unknown) angles PA ma j- n and 
all e $ yields the probability distribution Pj(e$ j, to) of measuring e$ in galaxy j, which is now only 
a function of to and the data. We proceed by defining the likelihood, C, for the parameter to to be 

C(e )= J2 In(^K^o)) -N, (16) 

j = l..N \ / 
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where the sum is taken over all N sample galaxies and e$ j is the measured value for the j 
galaxy. Once the best fit to is found, we estimate its variance by finding the parameter region whose 
likelihood is within unity of the best fit value. 

This analysis is illustrated in Figure 8, where we show a histogram of the measured e$ along 
with the expected distributions of measurements for to = 0.025, 0.045, and 0.075. Maximizing £(fo) 
yields the following results: (1) if we assume a spurious ellipticity due to spiral arms of 0.03 (see 
§3.2), which is added in quadrature to the measurement error, we find e = 0.045to °o% (2) if we 
do not include this source of error, then more of the observed ellipticity becomes attributable to 
the underlying disk ellipticity and to = 0.0551q 021 (3) if we assume that our estimate of cos(i)#/ 
from the linewidths is only an upper limit, and that in reality the galaxies are distributed uniformly 
between cos(i)#/ and cos(i) = 1, then less of the observed distortions is due to inclination effects 
and to increases to 0.0551q 021 anc ^ (4) if we assume the same conditions as in (1), but eliminate 
the two most strongly distorted galaxies (NGC 1302 and NGC 2718), we find e = 0.0301q ois- 
Therefore, we conclude that the galaxy disks in our sample have small but finite ellipticity. The 
characteristic ellipticity of the potential, within the disk plane, is to = 0.0451q 021 anc ^ this result 
is quite insensitive to the details of the statistical treatment. The potential ellipticity, e$ = 0.1, 
inferred for the Milky Way by Kuijken and Tremaine (1994) is somewhat higher than the mean 
found here, but still consistent with our distribution for e$ determined from external galaxies. 

It is worth stressing again that the biggest source of uncertainty in this analysis is the cross-talk 
between two-arm spirals and disk ellipticity. Even though we made an attempt to account for this 
through simulations, it is conceivable that spiral arm effects are bigger than we estimated. In this 
case, our to estimates still can be considered as a conservative upper limit on the triaxiality of the 
potential. 



4.2.3 Kinematic Implications 

These ellipticity estimates have kinematic consequences from which we can make predictions about 
the radial and azimuthal motions in disk galaxies. Elliptical streamlines arising in triaxial disks will 
result in nonzero radial velocity components for the local standard of rest (LSR) at all points along 
the orbit except at apocenter and pericenter. The LSR is defined to be co-moving with the net 
streaming of an ensemble of stars. Using Equation 8, the rms (denoted by (•••)) radial velocity of 
the LSR, averaged over all azimuthal angles is 

Using characteristic values for our sample, v c = 200 km s _1 and to = 0.045, we find that along a 
streamline the radial velocity with respect to the galaxy center varies from ~ —9km s _1 to ~ 9km s _1 
with (vr) = 6.4 km s _1 . 

This implies that the value inferred for the radial motion of the LSR in the solar neighborhood 
(14km s _1 , Blitz and Spergel 1990) is consistent with expectations based on the observed elliptical 
disks for a "typical" spiral galaxy. However, our data are equally consistent with the non-detection 
(vr = —1 ± 9km s _1 ) found by Kuijken and Tremaine (1994). The elliptical streamlines also lead 
to variations of the azimuthal velocity along the orbit (see Eq. 7). The fractional variance of this 
velocity is given by 

= US) 

As pointed out by Franx and de Zeeuw (1993), this variation contributes to the scatter in the 
TF relation — even if the true inclinations can be determined perfectly. Because the measured 
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H I linewidth is presumed to be proportional to 2vt a ngentiai{f) I sm it W1 H vary among identical 
galaxies seen at the same inclination but at different angles f. This effect leads to a scatter (in 
magnitudes) of at least 

for the TF relationship given by Pierce and Tully (1992). A characteristic ellipticity of to = 0.045 
would then produce a scatter of 0.15 magnitudes. However, an immediate comparison to the scatter 
in the TF relation from observed samples of galaxies is complex. The K'(2.2//m) photometry and 
H I linewidths presumably probe triaxiality at different radii. While our data probe typically 3 — 9 
kpc, the H I flux arises typically from radii larger by a factor of two. (S. Rao, priv. comm.). In 
addition, most TF observations are intended for distance estimates, therefore the samples are rarely 
statistically complete and may be biased against morphological or kinematic peculiarities. 

4.3 How Lopsided are Galaxy Disks? 

It is common lore that some galaxies appear lopsided, e.g. Sandage (1961) commented in the "Hubble 
Atlas" on the asymmetry of M101, NGC 1637, and others. Nonetheless, the study of Baldwin, 
Lynden-Bell, and Sancisi (1980) appears to be the only systematic study of such asymmetries in 
the literature. These authors studied lopsidedness in the H I distributions of galaxies, discussed 
excitation mechanisms, and estimated the longevity of such modes. The data presented in Figures 1 
and 3 show that lopsidedness in stellar disks is also common. Note that m = 1 amplitudes of i> 0.2 
are found in the outer parts of about one third of our sample. Regardless of the cause or longevity 
of such m = 1 distortions, they, like disk ellipticity, will lead to a radial velocity of the LSR in these 
galaxies. If we assume, in analogy to the m = 2 case, that to first order the closed streamlines at 
radius Ro have the shape 

R{f) = i? (l- ^cos^), (20) 
then the azimuthal and radial velocities have the form 

v v {if) = v c (l + cos if). (21) 



and 



v R {f) = -Vc^-sinf, (22) 

respectively. Applying the continuity equation (Eq. 9) to these orbits shows that m = 1 streamlines 
have the same fractional offset, e; op , as the isophotes. For an exponential brightness profile this 
offset is related to the Fourier amplitude by 



ei op =2A 1 -—. (23) 

Averaging over Eq. 22 and using Eq. 23, we now express the expected rms motion of the LSR at 
radius R in terms of the observables Ai and R exp , 

We rewrite the equation for (vr) using characteristic values in our sample 

ii 



(vr) = 7.4 km 



200 km 



0.11 



2 ^Rexp 



R 



(25) 
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Consequently, the expected coherent radial motions from lopsided and elliptical distortions are of the 
same magnitude and both are comparable to the value inferred for the solar neighborhood (e.g. Blitz 
and Spergel 1990). 

These dynamical inferences are valid regardless of the origin and the longevity of the m = 1 
distortions, provided the distortion persists longer than one orbital time. However, the question of 
the longevity of lopsidedness is both important and puzzling. As first discussed by Baldwin et al. 
(1980), if lopsidedness is due to elliptic orbits in an axisymmetric logarithmic potential, any initial 
azimuthal alignment of the orbit apocenters between the radii of R\ and i?2 would wind into a 
leading spiral at a rate of 

fi„„i ={y/2-l)vJ — — — ) . (26) 

\ dinner -Pouter J 

Using values characteristic for our sample, v c = 200 km s~ , R ou ter — 1.5i2 inner and Router — 8kpC, 
the "wind-up time," 2ir Q~j nd , is only about 10 9 yrs. This estimate of the lifetime in conjunction with 
our observation of significant lopsidedness in about 1/3 of our sample suggests that the generation 
of lopsidedness has a typical time-scale of 3xl0 9 yrs. Alternatively, we could suppose that disk stars 
merely respond to a potential distortion of the form $(i?, <p) = $(i?) + (e/ op /2) cosy. Then all closed 
orbits are aligned eccentric ellipses of the form given in Eq. (13) and there is no winding problem. 
However, as in the m = 2 case, the stellar orbits "oppose" the potential they reside in and cannot 
create such distortions self-consistently. Therefore the winding problem is merely replaced by the 
question of how to create a lopsided potential. 



4.4 How Strong are Spiral Arms? 

Two empirical statements made in Section 3.2 about the appearance of two-armed spirals at K' 
bear repeating : (1) they are usually coherent features, extending over a factor of about two in 
radius, and often winding through at least 180°, and (2) they are azimuthally broad features, in 
the sense that is the dominant even Fourier amplitude. Both of these properties are consistent 
with the hypothesis that the arms chiefly reflect variations in the stellar surface mass density (see 
Section 3.1). These IR observations are in support of optical studies (Schweizer, 1976; Jensen et 
a/.1981; Elmegreen and Elmegreen, 1984), which show insufficient color change across the arms 
to explain the intensity variations merely with stellar population changes. In Figure 9 we give a 
more quantitative picture of the spiral arm strength. In this Figure we plot the fractional variation 
A(i?) = Imax(R)/ Imin(R) — 1, where I m ax an d I m in are the maximum and minimum value of the 
surface brightness in a model image. To focus on features with two-fold symmetry we constructed 
a model image consisting only of the m = 0,2,4,6 Fourier components in the observed images. 
The uncertainties in Figure 9 were calculated from the uncertainties in the Fourier components. 
As discussed earlier, these intensity variations in azimuth reflect a mix of projection effects, disk 
ellipticity and spiral arms. However, projection effects contribute only (1 — cos i) R R and (disk) 

ellipticity accounts for only e$(l + R R v ) R R ; both effects cause contributions of A(i?) <^ 0.2. Since 
the variations displayed in Figure 1 are very often larger than these effects, they are most likely 
due to bars (at small radii) and spiral arms (at larger radii). Although we are not able to devise a 
physically well motivated definition of "spiral arm strength" , a comparison of Figure 9 with Figure 
2 (where one can check whether the m = 2 phase winds continuously) suggests that nearly half of 
the galaxies exhibit spiral arms with an arm-inter arm contrast, A(R), of about unity. 

If these luminosity variations indeed reflect mass variations, as the available evidence suggests, 
these spiral arms will have observable dynamical implications. We are able to make straightforward 
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predictions about non-circular motions for elliptical or lopsided disks because (as we showed in 
§4.2.1) the K' light traces the streamlines in a nonaxisymmetric potential. But even if we assume 
that the observed K' surface brightness //(i?, <p) perfectly traces the stellar surface mass density in 
the disk, T,(R, <p), and that the spiral pattern is stationary in some reference frame rotating at 
the analysis of the spiral arms is more complex than that of either ellipticity or lopsidedness (e.g. 
Contopoulos and Grosb0l, 1988, Grosb0l, 1993). 

Here we restrict ourselves to a brief discussion of the non-axisymmetric forces arising from the 
spiral arms, assuming that light traces mass. The Fourier decomposition of our images provides a 
convenient separation of the axisymmetric, Ao(R), and nonaxisymmetric, A m> o(R), components of 
the surface mass density. Using only Ao(R), we can determine the radial, f R (R), and azimuthal, 
/°, forces (which vanish identically) from the axisymmetric component of an observed galaxy. We 
calculate these forces by assuming a two-dimensional mass distribution and summing directly over 
the bins in our polar data grid. Similarly, we calculate ffi(R) and /™(i?) from n(R,<p). To focus 
on the effect of two armed spirals we include only the m = 2 and m = 4 components in our force 
calculation. The normalized forces arising from the nonaxisymmetric part, are defined by 

F R (R,V) = ^ (27) 

and 

^(ig,y)= ^y . (28) 

To account for centripetal force from an axisymmetric halo, the denominator should be replaced by 
f R (R) + ffi lo {R)- Here we calculate upper limits on Fr by adopting /^ a '°(i?) = 0. The nonax- 
isymmetric rms normalized forces, for m = 2 and 4, for two galaxies with strong spiral patterns, 
IC 2627 and NGC 6814, are between 5 and 15%. A similar analysis of I-band data for M81 (Visser 
1980), which is a weak-armed spiral (Schweizer 1976), produced an estimate of ~ 5%. This demon- 
strates that even though the mass density contrast across a spiral arm may be near unity, the force 
perturbations are a much smaller fraction of the total forces. Nevertheless, the forces are not in- 
significantly small and may pose problems for the standard linearity assumption in applications of 
the spiral density wave theory. 



5 Conclusions 

The study presented here shows that galaxy disks exhibit a bewildering variety of shapes even if 
one observes in the near-IR to minimize distortions due to dust and azimuthal mass-to-light ratio 
variations. We used Fourier techniques and model fitting to study the nonaxisymmetric component 
of the stellar disks. These techniques enabled us to decompose the nonaxisymmetric part of the 
K'(2.2//m) surface brightness into (a) disk ellipticity, (b) lopsidedness and (c) two arm spirals. The 
basic result from our analysis is that all three types of nonaxisymmetries are common. Our study 
indicates that distortions in the K'(2.2//m) light principally reflect distortions in the stellar surface 
mass density, rather than variations in the mass-to-light ratio. This mapping of light to mass can 
be tested further in dynamical studies or by comparing the current rate and distribution of star 
formation to the K'(2.2//m) light. 

Estimates of the dynamical effects arising from the various nonaxisymmetric components showed 
that all three types of distortions are comparably important, in the sense that they all induce veloc- 
ities (in stars or gas) which differ typically ~ 3 — 6% from those found in "equivalent" axisymmetric 



17 



galaxies. The peak values of the nonaxisymmetric motions are expected to be higher by factors of 
at least two and should be easily detectable in velocity maps of ionized or neutral gas. 

If we assume that the Milky Way is a "typical galaxy" , our results can provide a prior probability 
on local distortions. Even though we have not considered distortions with finite pattern speeds, it 
seems that all models suggested so far (Blitz and Spergel, 1990; Kuijken 1991; Kuijken and Tremaine, 
1994) have comparable prior likelihood. The one form of distortion that deserves more attention is 
massive spirals. Their modeling in the MW, however, is considerably more complex, especially since 
the spiral arms in other galaxies exhibit such a wide variety of morphologies. 

Our results confirm the assertion by Franx and de Zeeuw (1992) that non-axisymmetries may 
contribute significantly to the observed scatter in linewidth-luminosity relations for disk galaxies. 
Our results are most immediately applicable to studies where the kinematics are taken from ionized 
gas (e.g. Matthewson et a/.1992, Schommer et a/.1993), because similar radii are probed with 
these kinematics and our photometry. In these studies, disk ellipticities are expected to contribute 
~ 0.15mag scatter. 

The present results can be tested further in a variety of ways: the most obvious way is to check 
whether the non-circular motions suggested in §4 are borne out in the velocity fields of external 
galaxies. The most suitable kinematic data set would be HII velocity fields derived from Fabry- 
Perot imaging and HI 21-cm maps. For the first type of data it is often impossible to obtain a 
complete 2-D velocity filed, because the HII emission is very clumpy. For the second type of data it 
is difficult to obtain sufficiently high spatial resolution. Such a data model comparison also requires 
a more detailed prediction of the velocity field than provided in §4. For lopsided galaxies the most 
straightforward dynamical test is a search for a shift of the "systemic" velocity with radius. 

Finally, it will be interesting to explore the implications of the disk distortions for the shape 
of the dark halo. Throughout this paper we have only discussed the shape of the total potential. 
Since loop orbits "oppose" the shape of their potential, our results are consistent with halos that 
are quite distorted (in the disk plane); this is true in particular if the disks contribute much to the 
total gravitational force. 
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Appendix A. Model Fitting 



Here we describe in detail the parameterized model for the light distribution in the observed disk 
galaxies and the procedure used to find the best-fit parameters. This model is used in Section 3.2 
to determine the disk ellipticities. 

Al. The Bulge-Disk Model 

Our model for the light distribution of a sample galaxy consists of two components, a bulge and a 
disk, each of which is characterized by a parameterized fitting function. The bulge is assumed to 
have an i? 1 / 4 luminosity profile and constant ellipticity. The disk is assumed have an exponential 
luminosity profile and an ellipticity profile as expected for a constant potential ellipticity. 
The components are modeled with simple analytic expressions. The intensity at any point (R, (f>) is 
given by 

I(R^) = I B (R,<I>,P) + Id(R,<I>,p) + I.ky, (Al) 

where p denotes the vector of the model parameters, the origin of the coordinate system is assumed 
to be at the galaxy center, and Is and Ijj are the bulge and disk models respectively. The bulge 
model is given by 

/i? B (i?,«^ 1/4 



I B (RA) = h - exp | -7.67 

where 



V Re 



1 



(A2) 



R B (R, <j>) = R- ~ e B ) sin 2 (</>) + — !— cos 2 (</>), (A3) 
V 1 — eg 

I e is the effective intensity, R e is the effective radius, is the ellipticity, and is the minor axis 
position angle. The disk model is given by 

I D( R^) = I oe *J-M^l\, (A4) 
where Iq is the central disk surface brightness, R exp is the exponential scale radius, Hp is given by 

R D (R, <f>) = R-(l + ii, (R) ■ cos[2(^ - <f> D )]J (A5) 

and 

Qso(R) = e* • f 1 + t^-) • (A6) 

The decrease of the apparent eccentricity with increasing radius is expected for a disk embedded 
in a potential with constant distortion (Franx and de Zeeuw 1992). The sky level, I s ky, should be 
retained as a free variable, because the galaxies have a radial extent comparable to the field covered 
by the detector. 

A2. Determining the Best Parameters 
A2.1. Rebinning the Image 
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The "best fit model" is defined as the set of parameters, fbest, which minimizes 



„ . f \ Vdata(R, 9) J 

all data points 

With the algorithm employed here to find fbest, this sum must be evaluated many times (~ 1000). 
If this sum is taken over all pixels of a ~ 250 x 250 image it contains 62500 terms. A factor of ~ 100 
in computing time can be saved if the image is rebinned onto a (\og(R)-<f>) grid. Little information 
is lost in such rebinning because images of face-on galaxies have a geometry which is compatible 
with cylindrical coordinates and because the models have no features that cannot be captured by 
a ~ 30 x 24 (log(i?)-c/>) grid. Because pixels with R < 6 pixels (= 2.5") are affected by seeing, 
we restrict the polar grid to R > 6 pixels. For R ^> 10 pixels the (\og(R)-<f>) bins are statistically 
independent and \ 2 tests may be applied. A grid of 30 x 24 was found to have sufficient resolution 
ensure that neither the model nor the data vary significantly within a bin. 

A2.2. Biased Random Walk 

Because the parameter space spanned by the model has a high dimension and is infinite, the problem 
of finding the global minimum of \ 2 at pbest is difficult. Standard techniques of starting at an initial 
parameter guess, Pinit, and using a downhill gradient search (e.g. Press et al. 1986) are unreliable 
because it is unlikely that from a given point pi n u the path along the local \ 2 gradient will lead 
to fbest, rather than to a local extremum. Furthermore, we want to find the best solution for the 
parameters, subject to certain constraints, such as I e ,Io, R e , R exp ,eB,tD > 0. It is difficult to 
implement such constraints in gradient search methods. 

To avoid these difficulties, we use a biased random walk technique (often called simulated annealing, 
see Press et al. 1986), described in Rix and White (1992), that has proven to be a much more robust, 
though computationally expensive, means of finding the global minimum of \ 2 , independent of the 
starting parameter guess, Pinit- The application of the technique begins by selecting a "box" in 
parameter space defined by setting p™ n and pf' ax for all parameters on the basis of prejudice of 
permissible parameter ranges or physical constraints. Trial steps, bp, are then taken in parameter 
space starting from Pinit- At each step, \ 2 is evaluated and compared to \ 2 at the previous point 
in parameter space. The components of the steps, 8pi, are drawn randomly between —bpf' ax and 
bpf' ax . The vector of maximal steps, 8p max , must be specified. Such a step is always accepted 
in the approach to Xmin if ^X 2 < 0- However, the step is also accepted with a probability of 
exp(— Ax 2 /Axo) if > 0, where A%o is an adjustable "uphill penalty." All steps that would 

cause an exit of the parameter box are rejected. For a large initial value of \ 2 a few hundred to a 
few thousand steps are necessary. Once in the neighborhood of Xmin' the \ 2 penalty is stiffened and 
the maximal step sizes are decreased (see also Rix and White 1992). Finally, Xmin i s approached with 
a conventional gradient technique (Press et al. 1986) which enables us to calculate the co-variance 
matrix. 

Appendix B. The Age Distribution of Stars Observed in the B and K Bands 

Here we estimate the contribution of "young" stars to the total light of the galaxy at B and K(2.2;um) 
with a simple model. In this context, "young" means that t age , the age of the stars, satisfies 
tage ^ few tdynamical ~ 3 x 10 8 yrs. The K-band results are of immediate relevance to this paper, the 
B-band results are given for comparison. We use (Lq}\ to denote the mean luminosity at birth of a 
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set of stars for a given IMF, and s(t) to denote the rate at which stars, now at age t, were formed. 
At any given wavelength A, we define the fading of this set of stars with age as f\(t) such that 

L(t) x =f x (t){Lo)x, (CI) 

where L(t)\ is the luminosity at age t. The total luminosity due to stars younger than t can then 
be written as 

Lx(< t) = (L )x [ s(r)/ A (r)dr, (C2) 
Jo 

and the fraction of the luminosity they contribute to the total luminosity is 

R x (<t) = L x (<t)/L x (<t age ), (C3) 
where for a current galaxy t age = to ~ 10 10 years. 

To evaluate this integral exactly we need to know the star formation history and the fading function, 
which could be derived from stellar evolution models. To obtain a rough estimate of L\(< t), we 
make two simplifications: (1) we set s(t) = so + (2so/t age )t (Kennicutt 1983, Kennicutt et a/.1994); 
and (2) we set 

/o, fort<10 7 yrs; 



^ ~ \(</10 7 )-t, for*>10 7 yrs, 

where for the B-band /o ~ 1 and j ~ 0.84, for the K-band /o ~ 3 (to allow for the so-called 
red excursion of supergiant stars in the HR diagram), and j ~ 0.57. The latter quantities were 
extracted from the Figures given in Chariot and Bruzual (1991). We now express the luminosity of 
stars younger than t as 



L x (< t) = (L )x s 



t 1 -'* - 1 2(t 2 "T-l) 

Jo + 



(C4) 



1-7 (2 - j)t ag 

From this, Rx(< t) is easily calculated and is plotted in Figure 6 for both the B and K bands. 

Appendix C. Non-parametric Estimate of the Ellipticity Distribution 

Here we discuss how to estimate the most likely distribution of intrinsic disk ellipticities, P(e^), given 
a set of measured ellipticities, e$, with error Ae, for a galaxy observed at a kinematic inclination 
cos(i). The second angle, <f>, determining the disk orientation at fixed inclination i, is not known. 
Hence we must assume it to be distributed randomly. Let ej be the mean intrinsic ellipticity of the 
isophotes of the j th galaxy disk. Since most of the statistical weight in our fits of e$ j arises from 
radii R exp < R < 3R exp , we assume that ej = 1.5e$j (see Equation 11). An isophote can then be 
labeled by the parameter ip as 

r = cos(t/>)3? + (1 — ij) sin(i/>)j/ (Bl) 

and the viewing direction can be written as 

n = cos (<f>) sin(i)S? + sin(<^) sin(i)j/ + cos(i)z. (B2) 

Hence, at any point ip along the streamline the projected radius vector is given by r pro j = f x n 
with length 

r lroj = cos 2 (i) cos 2 (i/') + (l-?j) 2sm2 (V') +sin 2 (i) cos(V>) sin(</>) — (1 — ?j) sin(V>) cos(</>) . (53) 
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The ellipticity of the ellipse formed by the streamline is 



Ij = 1 — r pro j (min) / r pro j (max) . 



(54) 



We expect to measure (even without noise) e$j = fj/1.5. In the presence of a finite measurement 
error Aej, the distribution, 5(e$j), of expected measurements is given by 



This distribution is called the Rice distribution (see Wardle and Kronenberg 1974). For ?$ j Afj 
it tends to a Gaussian about ?$ j with dispersion Aej. When the underlying ellipticity is small 
compared to the measurement error, e$j <C Aej, the expectation value for a measurement does 
not tend to zero but rather to Afj . Note that the expression in B5 depends on the angles <f> 
and i through ?$j(<^,i). Averaging this distribution (B5) over all angles <f> yields the probability 
distribution, for the intrinsic ellipticity for the j th object. The final probability distribution, 

5(4), is obt ained by averaging Pj(e$) over our sample of galaxies, and this distribution is shown in 
Figure 7. Confidence regions for Pj(t$) were estimated by using the bootstrapping technique (Press 
et al. 1993). The technique consists of drawing 18 random members from the sample, allowing for 
multiple draws of individual sample members, and repeating the above analysis. By repeating many 
times this allows one to build an ensemble of possible p(e;„t r ) distributions with the same statistics as 
an ensemble drawn from the true p(e;„t r ). From the distribution of p(e 8ntr ), we calculate confidence 
limits on p(e;„t r ) for each value of e;„t r . The derived curve and error bounds are shown in Figure 7. 




(55) 
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TABLE 1 : The Sample 



Name 


Type 


mj 


D(') 

V / 


M B 


vti rfkm/sl 


WUkm/s] 

ZU L /J 


i(°) 

V / 


t e (s) 


NGC 600 


SBd 


12.9 


3.3 


-19.0 


1843 


80 


17 


1800 


NGC 991 


Sc 


12.4 


2.7 


-19.2 


1535 


82 


17 


1800 


NGC1015 a 


SBa 


13.0 


2.6 


-19.7 


2630 


87 


15 


1800 


NGC1302* 


Sa 


11.6 


4.8 


-20.2 


1703 


105 


16 


1620 


NGC1309* 


Sbc 


12.0 


2.2 


-20.3 


2135 


156 


23 


1800 


NGC1325A 


Sb 


13.3 


2.7 


-17.9 


1333 


46 


14 


1800 


NGC1376* 


Scd 


12.8 


2.0 


-21.0 


4162 


179 


21 


1800 


NGC1642* 


Sc 


13.3 


2.1 


-20.6 


4621 


143 


19 


1800 


NGC1703* 


SBb 


11.9 


4.5 


-19.6 


1526 


65 


11 


1980 


NGC2466* 


Sc 


13.5 


1.8 


-20.7 


5364 


175 


23 


1800 


NGC2485 C 


Sa 


13.1 


1.6 


-20.8 


4612 


200 


25 


1620 


NGC2718 


Sab 


12.7 


2.1 


-20.8 


3843 


130 


16 


1800 


NGC6814 


Sbc 


12.1 


3.0 


-19.5 


1565 


112 


21 


1080 


NGC7156 


Scd 


13.1 


1.6 


-20.5 


3984 


125 


17 


1500 


NGC7309 


Sbc 


13.0 


1.9 


-20.6 


4000 


136 


18 


1800 


NGC7742 


SO 


12.4 


2.0 


-19.4 


1649 


85 


17 


2160 


IC 2627 


Sbc 


12.6 


4.0 


-19.6 


2082 


42 


8 


360 


ESO-436 


Sc 


13.4 


1.8 


-20.2 


4079 


85 


13 


1260 



a HI profile remeasured by Tully (see Huchtmeier and Richter, 1989). 
' I-band image was obtained. 
c HI profile has broad wings. 
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TABLE 2 : Results from Fourier and Model Fitting Analysis 



Name 61 i? e:C p(arcsec) Ai(2.bR exp ) e$ PA(°) 



NGG 600 


17 


.6 


± 


2 


.8 





.11 


± 





.08 





.09 


± 





.06 


177 


± 


18 


NGG 991 


19 


.5 


-|- 


2 


.0 





.09 


-|- 





.05 





.10 


-|- 





.03 


89 


-|- 


09 


NGG1 01 5 


15 


2 


-|- 


1 


.5 





.23 


-|- 





.08 





.04 


-|- 





.03 


116 


-|- 


08 


NGG1 302 


11 


2 


-1- 





.8 





.03 


-1- 





.02 





.18 


-1- 





.03 


171 


-1- 


03 


NGG1 309 


11 


I 


-1- 





.6 





.37 


-1- 





.07 





.09 


-1- 





.01 


11 


-1- 


03 


NGG 1325 A 


21 


.5 


± 


3 


.0 





.58 


± 





.25 





.06 


± 





.07 


163 


± 


55 


NGG1376 


11 


.4 


± 





.6 





.05 


± 





.08 





.04 


± 





.03 


152 


± 


23 


NGG1 642 


9. 


4 


-1- 





\ 





.01 


-1- 





.02 





.04 


-1- 





.01 


117 


-1- 


08 


NGC1703 


18 


.0 


± 


1 


.5 





.17 


± 





.08 





.06 


± 





.04 


147 


± 


16 


NGC2466 


7. 


.7 


± 





.3 





.17 


± 





.04 





.10 


± 





.02 


7 


± 


02 


NGC2485 


14 


.5 


± 


2 


.4 





.27 


± 





.10 





.11 


± 





.05 


174 


± 


10 


NGC2718 


14 


.1 


± 


3 


.3 





.00 


± 





.07 





.27 


± 





.07 


156 


± 


06 


NGC6814 


18 


.6 


± 





.6 





.19 


± 





.04 





.09 


± 





.02 


9 


± 


05 


NGC7156 


9. 


.2 


± 





.4 





.08 


± 





.04 





.06 


± 





.02 


40 


± 


11 


NGC7309 


12 


.7 


± 





.4 





.19 


± 





.05 





.12 


± 





.02 


153 


± 


04 


NGC7742 


6. 


.6 


± 





.2 





.05 


± 





.01 





.02 


± 





.01 


162 


± 


15 


IC 2627 


21 


.2 


± 


1 


.2 





.17 


± 





.04 





.07 


± 





.03 


28 


± 


10 


ESO-436 


14 


.8 


± 


1 


.8 





.16 


± 





.14 





.10 


± 





.07 


164 


± 


20 
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Figure Captions 



Figure 1: K' images of the 18 galaxies in the sample. The angular size of each frame is 1.73' x 1.73'. 
The galaxies are shown in the same order as listed in Table 1 from left to right and top to bottom. 
North is at the top and East is to the left in all frames. Except for NGC 1302 and NGC 6814 all 
galaxies are displayed with the same stretch. 

Figure 2: Amplitudes and phases of the principal Fourier components in the K' images of the sample 
galaxies. Radius is plotted in disk scale-length units. Note that the amplitude scale on the bottom 
three panels has been expanded by 5/3. 

Figure 3: Histogram of disk lopsidedness for galaxies in the sample. The Fourier power Ai measured 
at 2.5 exponential disk scale lengths is expressed along the abscissa. The values of Ai have been 
corrected for their measurement errors (see text). 

Figure 4: Comparison between I and K' results for the six galaxies for which I-band data was 
obtained. The K' —I color, and the phases and amplitudes of the m = 1 and 2 components are 
shown (solid circles represent I data, open circles represent K' data). 

Figure 5: Model galaxy image, with m = 2 and 4 spirals characteristic of NGC 1703. The axisym- 
metric input disk and the arms result in an estimated disk ellipticity of only 0.03. 

Figure 6: The age distribution of the stellar K'(2.2 / um)-flux (open symbols) and B -flux (filled 
symbols), as predicted by the models described in Appendix C based on Chariot and Bruzual 
(1991). The top panel shows the predicted fraction of the total flux at current time {t age = 10 10 
yrs) arising from stars younger than t. Note that even in K'(2.2//m) the fraction of light from stars 
younger than a few dynamical periods, t ~ 3 x 10 8 yrs, may be 10%. The bottom panel shows the 
fraction of the luminosity from stars younger than t that arises from stars younger than 0.5t. 

Figure 7: The distribution of potential ellipticities P(e^). The solid line delineates the dis- 
tribution calculated from the data and the long dashed lines delineate the region of 95% confidence 
as estimated by "bootstrapping" (cf. Appendix B). The dotted line delineates the expected distri- 
bution from axisymmetric model disks at the inclinations given in Table 1 and measured with the 
observational errors estimated for our sample. 

Figure 8: The histogram of measured ellipticities, e$, and the expected distribution of such mea- 
surements for three assumed intrinsic ellipticity distributions are shown. The solid line assumes that 
P(e^,eo) = 1/eo exp(— e$/eo) with cq = 0.045. The more highly peaked dashed line illustrates the 
result from the same functional form but with to = 0.025, the other dashed line the result from a 
model with to = 0.075. 

Figure 9: Azimuthal intensity variations, as characterized by (T, max — S m8 ' n )/S m8 ' n . In order to 
focus on variations with two-fold symmetries, the maximal and minimal surface brightness were not 
calculated from the image itself, but from its Fourier components with m = 0, 2, 4, 6. The majority 
of galaxies show azimuthal variations in excess of unity at least at some radii. In some objects, e.g. 
NGC 1015, such strong variations are due to a bar, for many others, e.g. NGC 1309, NGC 2485, 
NGC 6814, NGC 7309, they are predominantly due to two-arm spirals. 
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